Understanding mobility in a social petri dish 
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Despite the recent availability of large data sets on human movements, a full understanding of 
the rules governing motion within social systems is still missing, due to incomplete information on 
the socio-economic factors and to often limited spatio-temporal resolutions. Here we study an entire 
society of individuals, the players of an online-game, with complete information on their movements 
in a network-shaped universe and on their social and economic interactions. Such a "socio-economic 
laboratory" allows to unveil the intricate interplay of spatial constraints, social and economic factors, 
and patterns of mobility. We find that the motion of individuals is not only constrained by physical 
distances, but also strongly shaped by the presence of socio-economic areas. These regions can be 
recovered perfectly by community detection methods solely based on the measured human dynamics. 
Moreover, we uncover that long-term memory in the time-order of visited locations is the essential 
ingredient for modeling the trajectories. 



Understanding the statistical patterns of human mo- 
bility, predicting trajectories and uncovering the mech- 
anisms behind human movements [l| is a considerable 
challenge with irnportant practical applications to traf- 
fic management 0, planning of urban spaces 0, Q, 
epidemics [BJsj , information spreading [lol [Tl| , and geo- 
marketing |12[. [isj . In the last years, advanced dig- 
ital technologies have provided huge amounts of data 
on human activities, allowing to extract information on 
human movements. For instance, observations of ban- 
knote circulation [3, mobile phone records fl6| . 
online location-based social networks [13, 18 1, GPS lo- 
cation data of vehicles [l^, or radio frequency identi- 
fication traces [l|, 0, [2^, have all been used as prox- 
ies for human movements. These studies have provided 
valuable insights into several aspects of human mobility, 
uncovering distinct features of human travel behaviour 
such as scaling laws IJ, |2l|, predictability of trajec- 
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22|, and impact of motion on disease spreading 
0-0, 23 1 . However, from a comparative analysis of the 
different works it emerges clearly that a "unified the- 
ory" of human mobility is still outstanding, since results, 
even on some very basic features of the motion, often 
appear to be contrasting One example is the mea- 
sured distribution of human trip lengths in various types 
of transportation: some studies agree that mobility is 
generally characterized by fat-tailed distributions of trip 
lengths [3, [2l| . while others report exponential or bi- 
nomial forms[l|, 0, [l^ . The discrepancies arise due to 
the different mobility data sets used, where mobility is 
indirectly inferred from some specific human activity in 
a particular context. For instance, mobile phone records 
typically provide location information only when a per- 




FIG. 1. The universe map of the massive multiplayer on- 
line game Pardus. The universe of Pardus can be represented 
as a network |24|1 with N = 400 nodes, called sectors (playing the 
role of cities), and K = 1160 links. Sectors are organized into 20 
different regions, called clusters, shown in the figure as different 
colour-shaded areas. There is no explicit set of goals in the game. 
Players are free to interact in a number of ways to e.g. increase their 
virtual wealth or status. Players move between sectors to interact 
with other players, e.g. to trade, attack, wage war, or to explore 
the virtual world. 



son uses the phone [21[ , while radio frequency identifica- 
tion traces like the ones of Oyster cards in the London 
subway Q only log movements based on public trans- 
portation systems. Analyses of these data sets can then 
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result in a possibly biased view of the underlying mobil- 
ity processes. Furthermore, most of the analyzed data 
sets have poor information on how socio-economic fac- 
tors influence human mobility patterns. More generally, 
the lack of an all-encompassing record set with positional 
raw data, including complete information on the socio- 
economic context and on the behaviour of all members 
of a human society, has so far limited the possibilities for 
a comprehensive exploration of human mobility. 

Here, we address the issue of mobility from a novel 
point of view by analyzing, with unprecedented preci- 
sion, the movements of a large number of individuals, 
the players of a self-developed massive multiplayer online 
game (MMOG). Such online platforms provide a fasci- 
nating new way of observing hundreds of thousands of 
interacting individuals who are simultaneously engaged 
in social and economic activities. The potential of on- 
line worlds as large-scale "socio-economic laboratories" 
has been demonstrated in a number of previous studies 
[ii-iil. For the MMOG at hand we have access to 
practically all actions (soj . including movements, accu- 
mulated over several years. This MMOG can therefore 
be considered as a "socio-economic petri dish" to study 
mobility in a completely controlled way. We can in fact 
observe the long-time evolution of a social system at the 
scale of an entire human society, having a perfect knowl- 
edge of all the spatio-temporal and socio-economic de- 
tails. In contrast to traditional studies in social science 
which are typically biased by well-known "interviewer- 
effects", in MMOGs the socio-economic measurements 
are objective and unobtrusive, since subjects are not con- 
sciously aware of being observed. 

Usi ng p ositional data of the players in the game uni- 
verse [3l| . in combination with other socio-economic in- 
formation from the game, we uncover various fundamen- 
tal features of mobility, and we provide a complete de- 
scription of the mechanisms causing the observed anoma- 
lous diffusion. Two are the main results of our work. 
First, we find the emergence of different spatial scales, 
due to the strong tendency of the players to limit their 
economic activities to some specific areas over long time 
periods and to avoid crossing the borders between dif- 
ferent areas. Making use of this observation, we propose 
an efficient method to identify socio-economic regions by 
means of community detection algorithms based solely 
on the measured movement dynamics. Our second re- 
sult unveils the driving mechanism behind the movement 
patterns of players: Locations are visited in a specific or- 
der, leading to strong long-term memory effects which 
are essential to understand and reproduce the observed 
trajectories. Finally, we provide large-scale evidence that 
neglecting either of these spatial or temporal constraints 
may obstruct the possibility of understanding the pro- 
cesses behind human mobility. 



A social arena: the online game Pardus 

Pardus is a massive multiplayer online game running 
since 2004, with a worldwide player base of more than 
350,000 individuals. It is an open-ended game whose 
players live in a virtual, futuristic universe and interact 
with each other in a multitude of ways. The topology 
of the universe can be represented as a network with 
400 nodes, called sectors, embedded in a two-dimensional 
space, the so-called universe map shown in Fig. [1] Each 
sector is like a city where players can have social rela- 
tions (establish new friendships, make enemies and wage 
wars), and entertain economic activities (trade and pro- 
duction of commodities). Typically, sectors adjacent on 
the universe map, as well as a few far-apart sectors, are 
interconnected by links which allow players to move from 
sector to sector. At any point in time, each sector is usu- 
ally attended by a large number of players. The network 
is sparse and, similarly to other spatial networks, is not a 
small world. It has a characteristic path length L ~ 11.89 
and a diameter dmax = 27, which means that, on average, 
players have to move through a non-negligible number of 
sectors to traverse the universe. See [Sj] Section S3 and 
js3 | Table I for a detailed characterization of the universe 
network structure. 

The sectors have been originally organized by the de- 
velopers of the game into 20 different clusters, which are 
perceived by the players as different political or socio- 
economic regions such as countries. For example, a player 
who is member of a political faction in the game is pro- 
vided some game-relevant protection in all clusters which 
are controlled by the faction, and has the opportunity 
of social promotion when accomplishing certain tasks 
within these clusters. Each cluster is shown in Fig. [T] 
with a different background colour. All clusters contain 
about 20 sectors each, with the exception of the central 
cluster, consisting of just one sector, and its surrounding 
three clusters having only 6-7 sectors. Sectors belong- 
ing to the same cluster are geographically close on the 
map, meaning that the distance between any two sectors 
in the same cluster is small, with an average distance 
around 3. Players typically have a "home cluster" where 
they focus their socio-economic activities over long time 
periods. Occasionally, they also move to sectors belong- 
ing to other clusters in order to explore the universe, to 
relocate their home (migrate), or during extreme game 
events such as wars. 

In Pardus, players are free to pursue whichever role 
they like to take. Many of them focus on expanding their 
social relations or political influence, some play the role 
of "scientists" exploring the universe, while others choose 
their main goal in trade and optimizing the amount of 
virtual money earned (25| . The large variety of complex 
socio-economic behaviours emerging in this online soci- 
ety, results in high heterogeneity in the mobility patterns. 
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in Figure [2] (a). For d < 15, the distribution is well-fitted 
by an exponential: 
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FIG. 2. Distribution of jump distances and of waiting 
times. To each player a time series consisting of the sector po- 
sitions over 1000 days is associated. A jump is said to occur when 
the sector position in the time series changes from one day to the 
following. The length d of a jump is measured in terms of graph 
distance and can take an integer value between 1 and dmax = 27, 
the diameter of the network, (a) The probability distribution of 
jump distances is reported in a semi-log plot. For d < 15, the dis- 
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tribution follows an exponential P{dj ~ e a with a characteristic 
length A Si 3. Players can also remain in the same sector for more 
days, without moving to other sectors. We define as waiting time 
At the number of consecutive days a player spends in only one sec- 
tor, (b) We show the probability distribution of waiting times At 
in a log-log plot, which is well fitted by a power-law P{At) ~ At"'' , 
with 13 « 2.2. 



such as observed in real human motion. However, differ- 
ently from other empirical studies on human movements, 
mobility in Pardus can be investigated in a controlled 
way, since complete information on actions of players is 
available [2^ [2a] ■ In this article we consider a data set 
consisting of movements in the network universe of all 
players who were active over a period of 1,000 days, as 
well as of socio-economic information about their envi- 
ronment. This opens the possibility of investigating mo- 
tion in relation to other social and economic factors. Note 
that we do not have to address the common issues of re- 
lying on incomplete data, on data that are only a proxy 
of mobility, or on data that ar e ag gregates of different 
types of transportation [§] . See [s3] Section 82 for more 
details on the data set. 



Basic features of the motion 

The position of each player in the universe, namely 
the ID number of the sector where the player is currently 
situated, is logged once a day. In this way the motion 
of each player becomes a time series of 1,000 sector po- 
sitions. A jump occurs when a player's sector position 
changes from one day to the following. The associated 
length d of a jump is measured in terms of graph dis- 
tance, an integer value between 1 and dmax = 27. The 
probability distribution of jump distances, computed for 
all players over the whole observation period, is reported 
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with a characteristic jump length A w 3. The existence of 
a typical travel distance, as also recently found in other 
mobility data [Hl,!!^, is related to the use of a single trans- 
portation mode in Pardus 32| . This allows to disentangle 
the intrinsic heterogeneity of the players from the effects 
due to the presence of different means of transportation 
0, which might be the cause of the scale- free distribu- 
tions found in mobile phone or other mobility data sets 



3, 16j. It has in fact been suggested that power laws 
in distance distributions of movement data may emerge 
from the coexistence of different scales [H, [H] . 

In some cases, players stay in the same sector for a 
number of consecutive days. For instance, 11 of the 1458 
considered players, although being active in the game, 
never jump within the entire observation period. On av- 
erage, a player does not change sector in approximately 
75% of the days. To better characterize the motion, we 
computed the waiting times At (measured in terms of 
number of days) between all pairs of consecutive jumps, 
over all players. The distribution of these waiting times, 
shown in Fig. [2] (b) follows a power-law distribution: 
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with an exponent 13 « 2.2, in agreement with other recent 
measurements on human dynamics (s^ . In addition, we 
found that the average waiting times of individual play- 
ers are distributed as a power-law (see [s3] Fig. 2). This 
implies a strong heterogeneity in the motion of differ- 
ent players, which is related to the heterogeneity in their 
general activity (see [sl] Section SI and [s3| Fig. 1). 



Mobility reveals socio-economic clusters 

Mobility patterns are influenced by the presence of 
the socio-economic regions in the network, highlighted 
in colours in Fig. [TJ The typical situation is illustrated 
in Fig. [3] (a) , with jumps within the same cluster being 
preferred to jumps between sectors in different clusters. 
In order to quantify this effect, we report in Fig. [3] (b), 
blue circles, the observed number of jumps of length d 
within the same cluster, divided by the total number of 
jumps of length d. This ratio is a decreasing function 
of the distance d, and reaches zero at d = 12, since no 
sectors at such distance do belong to the same cluster. 
As a null model we report the fraction of sector pairs 
at distance d which belong to the same cluster, see red 
squares in the same figure. The significant discrepancy 
between the two curves indicates that players indeed tend 
to avoid crossing the borders between clusters. For ex- 
ample, a jump of length d = 8 from one sector to another 
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clusters originally introduced by the developers of the 
game can be derived solely from the structure of the net- 
work. For this reason we adopted standard community 
detection methods based on the adjacency and on the 



distance matrix [36|, |37[. The results, reported respec- 



tively in Fig. m (d) and (e) , show that detected communi- 
ties deviate significantly from the clusters, implying that 
in our online world the socio-economic regions cannot 
be recovered merely from topological features. In com- 
parison we considered the player transition count matrix 
M, shown in Fig. 2] (c), which displays a similar block- 
diagonal structure as A and D, but with the qualitative 
difference that it contains dynamic information on the 
system. Figure |4] (f) shows that community detection 
methods applied to the transition count matrix M reveal 
almost perfectly all the socio-economic areas of the uni- 
verse. This finding demonstrates that mobility patterns 
contain fundamental information on the socio-economic 
constraints present in a social system. Therefore, a com- 
munity detection algorithm applied to raw mobility in- 
formation, as the one proposed here, is able to extract 
the underlying socio-economic features, which arc instead 
invisible to methods based solely on topology. For a de- 
tailed treatment of adopted community detection meth- 
ods and measures see [sl] Section S4, [sl] Table II and 
[sl Figs. 4 and 5. 

A long-term memory model 



FIG. 3. Influence of socio-economic clusters on mobility. 

(a) Sketch of jump patterns from a sector t to sectors within the 
same cluster, j and I, and to sectors in a different cluster, j' , V . Al- 
though sectors j' and V have the same graph distance from sector 
i as sectors j and I respectively, transitions across cluster border 
have smaller probabilities, (b) Quantitative evidence of the ten- 
dency of players to avoid crossing borders. Red squares show the 
null model, i.e. the fraction of all pairs of sectors at a given dis- 
tance d being in the same cluster. Blue circles show the fraction of 
measured jumps leading into the same cluster, per distance. Co- 
incidence of the two curves would indicate that clusters have no 
effect on mobility. Clearly this is not the case — there is a strong 
tendency of players to avoid crossing the borders between clusters. 



sector in the same cluster is expected only in 3% of the 
cases, while it is observed in about 20% of the cases. 

Now, the propensity of a player to spend long time 
periods within the same cluster might be simply related 
to the topology of the network, as in the case of ran- 
dom walkers whose motions are constrained on graphs 
with strong community structures [s^ . Nodes belonging 
to the same cluster are in fact either directly connected 
or are at short distance from one another. This prox- 
imity is reflected in the block-diagonal structure of the 
adjacency matrix A and of the distance matrix £), re- 
spectively shown in Fig. U] (a) and (b). We have there- 
fore checked whether the presence of the socio-economic 



In order to characterize the diffusion of players over 
the network, we have computed the mean square dis- 
placement (MSD) of their positions, cr^(t), as a function 
of time. Results reported in Fig. [5] (a) indicate that, for 
long times, the MSD increases as a power-law: 

a^{t)^t^ (3) 

with an exponent v « 0.26. This anomalous subdiffusive 
behaviour is not a simple effect of the topology of the 
Pardus universe. In fact, as shown in Fig. [5] (b), gray 
stars, the simulation of plain random walks on the same 
network produces a standard diffusion with an exponent 
V K, 1 nj) to t K, 100 days, and then a rapid saturation 
effect which is not present in the case of the human play- 
ers. 

Insights from the previous section suggest that the 
anomalous diffusion behaviour might be related to the 
tendency of players to avoid crossing borders. We have 
therefore considered a Markov model in which each 
walker moves from a current node i to a node j with 
a transition probability tt^ = my / rnn, where is 
the number of jumps between sector i and sector j, as 
expressed by the transition count matrix M of Fig. |4] (c) . 
The probabilities tt^ are the entries of the transition 
probability matrix 11, which contains all the information 
on the day-to-day movement of real players, such as the 
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FIG. 4. Extracting communities from network topology and from mobility patterns. (a) The adjacency matrix A of the 
universe network, (b) the matrix D of shortest path distances, and (c) the matrix M of transition counts of player jumps. Each of the 
three matrices contains 400 X 400 entries, whose values are colour-coded. Sector IDs are ordered by cluster, resulting in the block-diagonal 
form of the three matrices. We have used modularity-optimization algorithms to extract community structures from the information 
encoded in the three matrices. Different node colours represent the different communities found, while the 20 different colour-shaded 
areas indicate the predefined socio-economic clusters as in Fig. [T] The displayed Fowlkes and Mallows index G [0, 1] quantifies the 
overlap of the detected communities with the predefined clusters. The closer is to 1, the better the match, see [Sjl Section S4. (d) 
Although information contained in the adjacency matrix A allows to find 18 communities, a number close to the real number of clusters, 
the communities extracted do not correspond to the underlying colour-shades areas = 0.68). (e) Extracting communities from the 
distance matrix D only results in 6 different groups (J^ = 0.49). (f) The 23 communities detected using the transition count matrix M 
reproduce almost perfectly the real socio-economic clusters [J^ = 0.96), with only a few mismatched nodes detected as additional clusters. 
For more measures quantifying the match of communities, see [sj Table II. 



preference to move within clusters, the length distribu- 
tion of jumps, as well as the tendency to remain in the 
same sector. Despite this detailed amount of information 
used (the matrix IT has 160,000 elements), the Markov 
model fails to reproduce the asymptotic behaviour of the 
MSD, see magenta diamonds in Fig. [S] (b). Since the 
model considers only the position of the individual at its 
current time to determine its position at the following 
time, deviations from empirical data appear presum ably 
due to the presence of higher-order memory effects [38| . 
For this reason we have considered the recently proposed 
preferential return model (2l| which incorporates a strong 
memory feature. The model is based on a reinforcement 
mechanism which takes into account the propensity of 
individuals to return to locations they visited frequently 
before. This mechanism is able to reproduce the observed 
tendency of individuals to spend most of their time in a 



small number of locations, a tendency which is also preva- 
lent in the mobility behaviour of Pardus players (see [sj 
Fig. 3). However, the implementation of the preferential 
return model on the Pardus universe network is not able 
to capture the scaling patterns of the MSD, as shown in 
Fig. [5] (b). The reason is that in the model the proba- 
bility for an individual to move to a given location does 
not depend on the current location, nor on the order of 
previously visited locations. Instead, we observe that in 
reality individuals tend to return with higher probabil- 
ity to sectors they have visited recently and with lower 
probability to sectors visited a long time before. Conse- 
quently a sector that has been visited many times but 
with the most recent visit dating back one year has a 
lower probability to be visited again than a sector that 
has been visited just a few times but with the last visit 
dating back only one week. 
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FIG. 5. Diffusion scaling in empirical data and simulated 
models, (a) The mean square displacement (MSD) of the positions 
of players follows a power relation iT^(t) ~ with a subdiffusive 
exponent v ^ 0.26. The inset shows the average probability P<-o(t) 
for a player to return after t jumps to a sector previously visited. 
The curve follows a power law P<-j(r) ^ r~" with an exponent 
of a 1.3 and an exponential cutoff. We report, for comparison, 
(b) the MSD for various models of mobility. For random walk- 
ers and in the case of a Markov model with transition probability 
/J2i"^ij we observe an initial diffusion with an expo- 
nent V ai 1 and then a rapid saturation of cr^{t), due to the finite 
size of the network. A preferential return model also shows satu- 
ration and docs not fit the empirical observed scaling exponent i^. 
Conversely, a model with long-time memory (Time Order Memory) 
reproduces the exponent almost perfectly. Such a model makes use 
of the empirically observed P<-j(t) while the Markov model and 
the preferential return model over-emphasize preferences to loca- 
tions visited long ago and do not recreate the empirical curve well. 
Curves are shifted vertically for visual clarity. 



To highlight this mechanism we measured the return 
time distribution in the jump-time series (see Methods). 
In particular, we extracted the probability P<^(t) for an 
individual to return again (for the first time) to the cur- 
rently occupied sector after r jumps. As shown in the 
inset of Fig. [5] (a), we found that the return time distri- 
bution reads 



(4) 



re-visiting probability of recently explored locations. In 
this way we can capture the long-term scaling properties 
of movements. Exactly these asymptotic properties are 
fundamentally relevant for issues of epidemics spreading 
or traffic management. 

This "Time Order Memory" (TOM) model incorpo- 
rates a power-law distribution of first return times, to- 
gether with a power-law distribution of waiting times and 
an exponential distribution of jump distances, as those 
observed empirically in Fig.[5J We show below that these 
ingredients are sufficient to reproduce the subdiffusive 
behaviour reported in Fig. [5] (a). The model works as 
follows: an individual stands still in a given sector for 
a number of days drawn from the waiting time distribu- 
tion, Eq. Then, the individual jumps. There are two 
possibilities: (i) with a probability v she returns to an 
already visited sector, (ii) with the probability 1 — v she 
jumps to a so far unexplored sector. In case (i), one of the 
previously visited sectors is chosen according to Eq. (jH). 
In the exploration case (ii), the individual draws a dis- 
tance d from the distance distribution, Eq. ([1]), and jumps 
to a randomly selected, unexplored sector at that dis- 
tance. The model has four parameters. The parameters 
A, P and a of equations (HJ, ([2]) and (|4]) respectively, are 
fixed by the data. Further, averaging over all jumps and 
players, the probability of returning to an already visited 
location is v w 0.83. Similarly to the measured data, 
the MSD of the TOM model, black squares in Fig.[5](b), 
exhibits no saturation effects and displays an exponent 
i^TOM = 0.23 ± 0.02 (the error is calculated over an en- 
semble of realizations) in agreement with the exponent 
observed for the players. 

The flat slope oi u ^ 0.26 and the lack of saturation 
of the MSD of the players over the whole observation pe- 
riod exposes the significant level of subdiffusivity in the 
motions of individuals, consistent with previous findings 
[2ll . [39l - l42| . However, the mere tendency of individuals to 
return to already visited locations is not sufficient to cap- 
ture these subdiffusive properties of the MSD, but it is 
fundamental to consider a mechanism that takes into ac- 
count the temporal order of visited locations, as achieved 
by the TOM model. Moreover, the TOM model is real- 
istic in the sense that, in contrast to Markov models, it 
takes into account the tendency of individuals to develop 
a preference for visiting certain locations. At the same 
time it allows for the possibility that a previously pre- 
ferred location becomes not frequented anymore. This 
view provides an alternative to recently suggested rein- 



forcement mechanisms in preferential return models [21 



with an exponent a « 1.3. We used this information for 
constructing a model which takes into account the higher 



The possibility for individuals to "change home" is rele- 
vant when the model should be able to account for mi- 
gration, which is an important feature in the long-time 
mobility behaviour of humans. 

Finally, we discuss to which extent the findings from 
our "social petri dish" are valid also for human popula- 
tions unrelated to the game. Previous analyses of human 
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social behaviour in Pardus [25|,|26[ have shown agreement 
with well-known sociological theories and with properties 
on comparable behavioural data. Examining the prefer- 
ence of players to move within socio-economic regions 
is of obvious importance for clearing up the role of po- 
litical or socio-economic borders on the movement and 
migration of humans, where the presence of borders has 
a strong influence on mobility [l5, 43-4^. Online soci- 
eties as the one of Pardus have the evident potential to 
serve as "socio-economic laboratories", where the com- 
plete knowledge of activities, social relations, and posi- 
tions of all individuals can significantly advance our un- 
derstanding of large-scale human behaviour, in particular 
of mobility. 



Methods 

Data set. We focus on one of the three Pardus universes, 
Artemis. For this universe, we extract player mobility data from 
day 200 to day 1200 of its existence. We discard the first 200 days 
because social networks between players of Pardus have shown ag- 
ing effects in the beginning of the universe, i.e. there seems to exist 
a transient phase in the development of the society, possibly affect- 
ing mobility, which we would like to avoid considering [25j . To 
make sure we only consider active players, we select all who exist 
in the game between the days 200 and 1200, yielding 1458 players 
active over a time-period of 1000 days. The sector IDs of these 
players, i.e. their positions on the universe network's nodes, are 
logged every day at 05:35 GMT. Players typically log in once a day 
and perform all their limited movements of the day within a few 
minutes, see [sJl Section SI. The legal department of the Medical 
University of Vienna has attested the innocuousness of the used 
anonymized data. 

Transition count matrix and transition probability ma- 
trix. The entry rriij of the transition count matrix M is equal to 
the number of times a player's position was on sector i and then, 
on the following day, on sector j. This number is cumulated for 
all players. The entry Tr^j of the transition probability matrix 11 
corresponds to the probability that a player moves to a sector j 
given that on the previous day the player's location was sector i. It 



movements from sector i to sector j, and the sum over / is over all 
sectors of the universe. The matrix 11 is a stochastic matrix, i.e. it 
has the property that the entries of each row sum to one. 

MSD and diffusion. The MSD is defined as cr^ {t) = 
((r (T + t) - r (T))^), where r (T) and r(T + t) are the sectors a 
player occupies at times T and T + t respectively, and where by 
{r {T + t) — r (T)) we denote the distance between the two sectors. 
The average (■) is performed over all windows of size t, with their 
left boundaries going from T=0 to T=1000-t, and over all the 1458 
players in the data set. If has the form <r^{t) ~ t" with an 
exponent u < I, the diffusion process is subdiffusive, in the case 
1/ > 1 it is super-diffusive. An exponent of ^ = 1 corresponds to 
classical brownian motion [sol . |40| . 

Jump-time and first return time distribution. We trans- 
form the time-series of daily sector IDs occupied by the play- 
ers from real-time to jump-time, in order to be able to com- 
pare time-series of different length and to focus on the move- 
ments between sectors. An example of this conversion is provided: 
a time series [5, 5, 5, 32, 32, 104, 5, 5, 104, 104, 104, 32, 337, 337, 32. . .] 
becomes in jump-time [5,32,104,5,104,32,337,32,...]. We de- 
note jump-time by the greek letter r, that is, at jump-time t a 
player has performed exactly t jumps. We use T in the compu- 
tation of the first return time distribution. In the hypothetical 
time series of sectors [5, 32, 104, 5, 104, 32, 337, 32] a first return to 



a sector lying r = 1 jumps back happens 2 times (104, 5, 104 and 
32,337,32), for t = 2 this happens once (5,32,104,5), for r = 3 
also once (32, 104, 5, 104, 32). Hence, in this example, P^{1) = 0.5, 
P<_j(2) = P<-j(3) = 0.25, where the sum over all P<-j(t) is equal to 
1. 



Acknowledgments 



This work was conducted under the HPC-EUROPA2 project 
(project number: 228398) with the support of the European Com- 
mission — Capacities Area - Research Infrastructures initiative, 
and within the framework of European Cooperation in Science 
and Technology Action MP0801 Physics of Competition and Con- 
flicts. M.S. and S.T. acknowledge support from the Austrian Sci- 
ence Fund Fonds zur Forderung der wissenschaftlichen Forschung 
P 23378, and from project EU FP7 - INSITE. M.S., R.S. and G.P. 
also thank the Santa Fc Institute for the opportunities offered dur- 
ing the Complex Systems Summer School 2010, where this project 
originated. 



[1] Barthelemy, M. Spatial networks. Phys. Rep. 499, 1-101 

(2010) . 

[2] Guimera, R., Mossa, S., Turtschi, A. & Amaral, L. 

The worldwide air transportation network: anomalous 

centrality, community structure, and cities global roles. 

Proc. Natl. Acad. Sa. USA 102, 7794-7799 (2005). 
[3] Helbing, D. Traffic and related self-driven many-particle 

systems. Rev. of Mod. Phys. 73, 1067 (2001). 
[4] Makse, H. A., Havlin, S. & Stanley, H. E. Modelling 

urban growth patterns. Nature 377, 608-612 (1995). 
[5] Roth, C, Kang, S. M., Batty, M. & Barthelemy, M. 

Structure of urban movements: Polycentric activity and 

entangled hierarchical flows. PLoS ONE 6, el5923 

(2011) . 

[6] Pastor-Satorras, R. & Vespignani, A. Epidemic spreading 
in scale-free networks. Phys. Rev. Lett. 86, 3200-3203 
(2001). 

[7] Colizza, v., Barrat, A., Barthelemy, M. & Vespignani, 

A. The role of the airline transportation network in the 

prediction and predictability of global epidemics. Proc. 

Natl. Acad. Set. USA 103, 2015 (2006). 
[8] Hufnagel, L., Brockmann, D. & Geisel, T. Forecast and 

control of epidemics in a globalized worlds. Proc. Natl. 

Acad. Sci. USA 101, 15124-15129 (2004). 
[9] Balcan, D. et al. Multiscale mobility networks and the 

spatial spreading of infectious diseases. Proc. Natl. Acad. 

Sci. USA 106, 21484-21489 (2009). 
[10] Miritello, G., Moro, E. & Lara, R. Dynamical strength 

of social ties in information spreading. Phys. Rev. E 83, 

045102 (2011). 

[11] Onnela, J. et al. Structure and tie strengths in mobile 
communication networks. Proc. Natl. Acad. Sci. USA 
104, 7332 (2007). 

[12] Querela, D., Lathia, N., Calabrese, P., Di Lorenzo, G. 
& Crowcroft, J. Recommending social events from mo- 
bile phone location data. In Data Mining (ICDM), 2010 
IEEE 10th International Conference on, 971 -976 (2010). 

[13] Jensen, P. Network-based predictions of retail store com- 
mercial categories and optimal locations. Phys. Rev. E 
74, 035101 (2006). 



8 



Brockmann, D., Hufnagel, L. & Geisel, T. The scaling 
laws of human travel. Nature 439, 462-465 (2006). 
Thiemann, C, Theis, F., Grady, D., Brune, R. & 
Dirk Brockmann, D. The structure of borders in a small 
world. PLoS one 5, el5422 (2010). 

Gonzalez, M., Hidalgo, C. & Barabasi, A. Understanding 
individual human mobility patterns. Nature 453, 779- 
782 (2008). 

Scellato, S., Noulas, A., Lambiotte, R. & Mascolo, C. 
Socio-spatial properties of online location-based social 
networks. Proceedings of ICWSM 11 (2011). 
Scellato, S., Musolesi, M., Mascolo, C., Latora, V. & 
Campbell, A. Nextplace: A spatio-temporal prediction 
framework for pervasive systems. Pervasive Computing 
152-169 (2011). 

Bazzani, A., Giorgini, B., Rambaldi, S., Gallotti, R. & 
Giovannini, L. Statistical laws in urban mobility from 
microscopic gps data in the area of florence. J. Stat. 
Mech. 2010, P05001 (2010). 

Cattuto, C. et al. Dynamics of person-to-person interac- 
tions from distributed rfid sensor networks. PloS one 5, 
ell596 (2010). 

Song, C., Koren, T., Wang, P. & Barabasi, A. Modelling 
the scaling properties of human mobility. Nature Physics 
6, 818-823 (2010). 

Song, C., Qu, Z., Blumm, N. & Barabasi, A. Limits 
of predictability in human mobility. Science 327, 1018 
(2010). 

Belik, v., Geisel, T. & Brockmann, D. Natural human 
mobility patterns and spatial spread of infectious dis- 
eases. Phys. Rev. X 1, 011001 (2011). 
Boccaletti, S., Latora, V., Moreno, Y., Chavez, M. & 
Hwang, D. Complex networks: structure and dynamics. 
Phys. Rep. 424, 175-308 (2006). 

Szell, M. & Thurner, S. Measuring social dynamics in 
a massive multiplayer online game. Social Networks 32, 
313-329 (2010). 

Szell, M., Lambiotte, R. & Thurner, S. Multirelational 
organization of large-scale social networks in an online 
world. Proc. Natl. Acad. Sci. USA 107, 13636-13641 
(2010). 

Castronova, E. On the research value of large games. 

Games and Culture 1, 163-186 (2006). 

Bainbridge, W. The scientific research potential of virtual 

worlds. Science 317, 472 (2007). 

www.pardus.at. 

Thurner, S., Szell, M. & Sinatra, R. Emergence of good 
conduct, scaling and zipf laws in human behavioral se- 



quences in an online world. arXiv:1107.0392vl (2011). 

[31] Whenever we address the position or the movement of 
'a player' in the game universe, this is meant as a short 
form for referring to the virtual avatar which is uniquely 
associated to and controlled by the player. This abbrevia- 
tion is consistent with the tendency of players to identify 
themselves with their avatars. 

[32] Kolbl, R. & Helbing, D. Energy laws in human travel 
behaviour. New J. of Phys. 5, 48 (2003). 

[33] Han, X., Hao, Q., Wang, B. & Zhou, T. Origin of the 
scaling law in human mobility: Hierarchy of traffic sys- 
tems. Phys. Rev. E 83, 036117 (2011). 

[34] Barabasi, A. The origin of bursts and heavy tails in 
humans dynamics. Nature 435, 207 (2005). 

[35] Fortunate, S. Community detection in graphs. Phys. 
Rep. 486, 75-174 (2010). 

[36] Arenas, A., Fernandez, A. & Gomez, S. Analysis of the 
structure of complex networks at different resolution lev- 
els. New J. of Phys. 10, 053039 (2008). 

[37] Newman, M. Analysis of weighted networks. Phys. Rev. 
E 70, 056131 (2004). 

[38] Sinatra, R., Condorelli, D. & Latora, V. Networks of 
motifs from sequences of symbols. Phys. Rev. Lett. 105, 
178702 (2010). 

[39] West, B., Grigolini, P., Metzler, R. & Nonnenmacher, T. 

Fractional diffusion and levy stable processes. Phys. Rev. 

E 55, 99 (1997). 
[40] Metzler, R. & Klafter, J. The random walk's guide to 

anomalous diffusion: a fractional dynamics approach. 

Phys. Rep. 339, 1-77 (2000). 
[41] Scafetta, N., Latora, V. & Grigolini, P. Levy statistics in 

coding and non-coding nucleotide sequences. Phys. Lett. 

A 299, 565-570 (2002). 
[42] Viswanathan, G. et al. Optimizing the success of random 

searches. Nature 401, 911-914 (1999). 
[43] Ratti, C. et al. Redrawing the map of Great Britain from 

a network of human interactions. PLoS One 5, el4248 

(2010). 

[44] Newman, D. The lines that continue to separate us: bor- 
ders in our borderless world. Progress in Human Geog- 
raphy 30, 143 (2006). 

[45] Lambiotte, R. et al. Geographical dispersal of mobile 
communication networks. Physica A 387, 5317-5325 
(2008). 

[SI] Supporting Information: |http : //www, complex- systems .meduniwien 



